Mon. Not. R. Astron. Soc. 000, 000-000 (1997) Printed 1 February 2008 (MN plain T^ macros vl.6) 

Magnetohydrodynamics in Superconducting-Superfluid 
Neutron Stars 



Gregory Mendell 

Institute for Fundamental Theory, University of Florida, Gainesville, FL 32611, U.S.A. 
Email: mendell@phys.ufi.edu 



a^ 

CD 

a^ 

m 
> 

m 

o 

(N 

O 



CL| 1 INTRODUCTION 



ABSTRACT 

MHD equations are presented for the mixture of superfluid neutrons, superconducting 
protons, and normal electrons believed to exist in the outer cores of neutron stars. 
The dissipative effects of electron viscosity and mutual friction due to electron- vortex 
scattering are also included. It is shown that Alfven waves are replaced by cyclotron- 
vortex waves that have not been previously derived from MHD theory. The cyclotron- 
vortex waves are analogous to Alfven waves with the tension due to the magnetic 
energy density replaced by the vortex energy density. The equations are then put 
into a simplified form useful for studing the effect of the interior magnetic field on 
the dynamics. Of particular interest is the crust-core coupling time which can be 
inferred from pulsar glitch observations. The hypothesis that cyclotron-vortex waves 
play a significant role in the core spin-up during a glitch is used to place limits on the 
interior magnetic field. The results are compared with those of other studies. 
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Below a superfluid transition temperature of approximately 
10^ K the outer core of a neutron star is predicted to exist as 
an exotic plasma of mostly superfluid neutrons, with a small 
fraction of type II superconducting protons, normal elec- 
trons, and even fewer muons (Baym, Pethick & Pines 1969; 
Baym, Bethe & Pethick 1971; Baym & Pethick 1975, 1979; 
Epstein 1988) . The mass density of this region ranges from 
2.8 X 10^* g • cm~'^ to a few times this value, and the protons 
and electrons make up only a few percent of the total parti- 
cle density. Rotation causes the superfluid neutrons to form 
an array of quantized vortices and magnetic fields cause the 
superconducting protons to form an array of quantized flux 
tubes (referred to generically as vortices from here on). Pul- 
sar observations show that neutron stars can rotate at least 
642 times per second and typically have surface magnetic 
fields of about 10^^ G (Lyne 1995; Nice 1995). Weaker fields 
in the range 10*- 10^° G have also been observed in neutron 
stars, mostly in binary systems. The evolution of the surface 
field and its link to the interior field can provide clues to the 
physics of the neutron-star core (Bhattachary 1995). Other 
clues come from thermal-evolution studies (Tsuruta 1995) 
and the observation of rapid spin-ups, or glitches, of several 
pulsars (Alpar & Pines 1985; Sauls 1989; Alpar 1995; Chau 
1995; Link & Epstein 1996.) 

In the past it has been suggested that magnetic fields 
are responsible for the core spin-up during a pulsar glitch 
and that glitches may excite Alfven waves in this region 
(Baym, Pethick, Pines & Ruderman 1969; Van Horn 1980 
and references therein). More recently, a constraint on the 
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crust-core coupling time for the Vela pulsar has been used 
to set lower limits on the magnitude of the interior mag- 
netic field (Abney, Epstein & Olinto 1996). These lower lim- 
its were derived using the Alfven wave velocity, though a 
'modified' Alfven wave velocity was used in the case of a 
superconducting-superfiuid interior . 

The first goal of this paper is to develop the mag- 
netohydrodynamic (MHD) theory of the superconducting- 
superfiuid outer-core of a neutron star, and to study the 
characteristic waves that occur in this exotic plasma. Hydro- 
dynamic equations for this region have been developed by 
MendeU & Lindblom (1991), MendeU (1991a,b), and Lind- 
blom & Mendell (1994). Here the MHD equations given in 
Mendell (1991a) are developed with a focus on the effects 
of the proton vortices. It is argued that the energy density 
due to individual proton vortices dominates the vortex ef- 
fects, and that the interactions between vortices, including 
pinning between the neutron and proton vortices, can be ig- 
nored. After the MHD equations have been presented it is 
shown that Alfven waves do not occur in the plasma consid- 
ered here, but are replaced by cyclotron-vortex waves that 
have not been previously derived from hydrodynamic the- 
ory. The cyclotron-vortex waves found here are shown to be 
equivalent to the 'modified' Alfven waves used in Abney et 
al. (1996). 

The second goal of this paper is to put the equations 
into a simplified form useful for studying the effects of the 
interior magnetic field on the dynamics. The main simplifica- 
tion is to drop all but the largest of the force terms due to the 
magnetic flux and the vortices. The resulting equations are 
the superconducting-superfluid generalization of the MHD 
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equations that Easson & Pethick (1979) gave for the case 
of an ordinary p-e plasma. One advantage of the simphfied 
equations is that the dissipative effects of electron viscosity 
and mutual friction due to electron- vortex scattering are eas- 
ily studied following the methods given in Mendell (1991b). 
The scattering of electrons off the neutron vortices usually 
dominates the dissipative effects, though electron viscosity 
can be important in some situations. Another advantage of 
the simplified equations is that they should be particularly 
convenient for studying the effects of the magnetic field and 
cyclotron-vortex waves on the core spin-up during a pulsar 
glitch. 

The final goal of this paper then is to investigate the hy- 
pothesis that cyclotron- vortex waves play a significant role 
in the crust-core coupling time, which can be inferred from 
pulsar glitch observations. Simple estimates are used to place 
several limits on the interior magnetic field, and these are 
compared with the results of other studies. 

The paper is organized as follows. In Section 2 the rel- 
evant properties of superconductivity and superfluidity are 
reviewed and the MHD equations for the outer-core region 
are presented. The characteristic frequencies associated with 
these equations are also reviewed. In Section 3 it is shown 
that Alfven waves do not exist and the dispersion relation for 
transverse axial modes in a constant magnetic field is found. 
The frequency of these modes depends on the cyclotron fre- 
quency and Kelvin vortex oscillation frequency, and so is 
called the cyclotron-vortex frequency. In Section 4 the ef- 
fects of rotation are included and it is shown that the usual 
inertial mode exists, along with an inertial mode modified 
by Fermi liquid 'drag' effects, and a mode dependent on the 
cyclotron-vortex frequency. In Section 5 the equations are 
simplified, and dissipative effects are explored. In Section 6 
the spin-up of the core during a pulsar glitch is discussed. In 
Section 7 the results are briefly summarized and suggestions 
for further work are discussed. 



to the MHD equations. A more complete review is given in 
Appendix A of Mendell (1991a). 

In all the equations that follow the subscripts n, p, and 
e refer to the neutrons, protons, and electrons respectively. 
For each species, p is the mass density, v is the velocity field, 
m is the mass, and e, when not a subscript, is the absolute 
value of the charge of the electron. Standard symbols are 
used for the fundamental constants fi, c, fee, and G, though 
G is also used as an abbreviation for Gauss, and Gaussian 
units are used throughout. 

As in a laboratory type II superconductor (e.g., see 
Tilley & Tilley 1986), each proton vortex will carry one 
quantum of flux 



$o = {hc/2e) 



(1) 



and each vortex has a core region of normal protons with a 
radius equal to the coherence length, given by 



c = 



Ahvi 



7-KkBT, 



^^ ~ 4.9 X 10-^M ^ ) Pr>.^"T, 



cp9 



(2) 



where v^p is the Fermi velocity, ppia is the mass density of 
the protons divided by lO^^g • cm^^, and Tcpg is superfluid 
transition temperature divided by 10^ K. Outside the vortex 
core the proton velocity field and magnetic induction fall off 
exponetially with a decay length-scale known as the London 
depth, given by 



A = 



IpC'^ 



4-Ke^Pp 



1.3 X 10" 



Pp 
Ppp 



1/2 



Ppi3 ~^'^ cm. 



(3) 



The assumption of type II superconductivity is valid when 
\/2A > ^. In these equations rrip is the effective mass of the 
protons, while ppp is a component of the mass density matrix 
[see the discussion around eqs. (11) and (12)]. The energy 
per unit length of a proton vortex, Sp, is then specified by the 
kinetic energy of the protons plus the energy of the magnetic 
induction, i.e., 



2 VORTICES AND THE MHD EQUATIONS 

An ordinary plasma is a mixture of oppositely charged flu- 
ids that is locally charge neutral when in equilibrium. High 
frequency plasma waves result when the system is perturbed 
away from charge neutrality. But low frequency waves are 
also possible and these waves can be studied by taking the 
MHD limit. This limit is in force when the following restric- 
tions apply: 1) the perturbations away from equilibrium pre- 
serve charge neutrality and, 2) the phase velocities are small 
compared to the speed of light. Alfven waves are the charac- 
teristic MHD waves that travel along magnetic field lines in 
an ordinary plasma. (For more on the MHD limit see Boyd 
& Sanderson 1969; Jackson 1975; Easson & Pethick 1979; 
Freidberg 1982.) 

Now consider the plasma in the outer-core region of 
a rotating and magnetized neutron star consisting of su- 
perconducting protons and normal electrons, immersed in a 
background of superfluid neutrons. A major difference be- 
tween this plasma and an ordinary one is that the neutrons 
and protons are expected to form arrays of quantized vor- 
tices (see especially Baym, Pethick & Pines 1969 concerning 
the proton vortices). Since the effects of the proton vortices 
play a dominate role in the results to be discussed here, it 
will be useful to review their properties before proceeding 



"Ppp^P 



Stt 



da. 



(4) 



where Vp is the velocity field of the protons, B is the mag- 
netic induction of the vortex, and the integral is understood 
to extend over the plane perpendicular to the axis of the 
vortex. As is shown in standard texts on superconductivity 
(and Mendell 1991a) Sp can be approximated as. 



<i>^ 



167r2A2 



ln(A/C) 



^2.1 X 10 
where 
A = 1.6 In 



^ ) (1 + A)ppi3 erg • cm-^ 

Pp 



(5) 



1/2 



-^ cp^PplS 



-5/6 



(6) 



This approximation is valid as long as A > ^ which implies 
A > —1; it is assumed here that A is small or at most of 
order unity. 

In a typical neutron star the spacing between vortices is 
much larger than the London depth but still much smaller 
than the radius of a neutron star. These statements apply to 
the neutron vortices as well as the proton vortices. Assuming 
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a triangular array in equilibrium (as is observed in terrestrial 
superfluids) the intervortex spacings between neutron and 
proton vortices, d„ and dp, are given respectively by 



rfn = 



rfn = 



h 



h 



1/2 



1/2 



-1/2 



3.5 X 10"'''Boi2"'^^ 



(7) 



(8) 



where 0,2 is the rotational angular velocity divided by 
10^ s^^, Boi2 is the average equilibrium magnetic induc- 
tion divided by 10^^ G, and Qc is the cyclotron frequency 
of the protons. 

To investigate the hydrodynamic behavior of the fluids 
on length-scales much larger than the intervortex spacing 
it is therefore appropriate to perform a smooth averaging 
over the vortices, as discussed e.g., by Baym & Chandler 
(1983) (for a review see Sonin 1987). Following these refer- 
ences the number of neutron and proton vortices per unit 
area, n^n and rivp, are related to the smooth-averaged vor- 
ticities and magnetic induction by the following form of the 
London equations: 



h 



2fn„ 



a„B = 



2mD 



(9) 



(10) 



In these equations Up = e/nipC, and the unit vectors i/„ and 
Vp point in the smooth-averaged direction of the vortices. 
In equation (10) B is the smooth-averaged magnetic induc- 
tion. The smooth-averaged vorticities in these equations are 
given by u)„ = V x u„ and cjp = V x tip, where t)„ and 
Vp are the smooth-averaged superfluid velocities. The Lon- 
don equations then can be taken as the definition of these 
velocities. Given these definitions it must be noted that due 
to Fermi liquid 'drag' effects the mass currents of the neu- 
trons and protons depend on the superfluid velocities of both 
species via the mass density matrix p„„, p„p — pp„, and 
ppp, as discussed by Andreev & Bashkin (1976) and Alpar, 
Langer & Sauls (1984). The recent calculation of Borumand, 
Joynt & Kluzniak (1996) shows that 



Pnp Ppn Pp 



1-^(1 + Ff73) 



pp - ppp 



(11) 



(12) 



where rrip and Ff^ are the effective mass of the proton and 
a generalized version of the Landau parameter respectively. 
These factors are defined in Fermi liquid theory in terms 
of quasiparticle interactions that must be determined either 
from nuclear matter theory or by observation. Most of the 
previous studies of superfluidity in neutron stars have ig- 
nored the factor Ff^. However, because of the uncertainties 
in rrip and Ff'^, all studies agree that for the neutron star 
outer core p„p ~ —pp, pnn ~ pn, and ppp ~ 2pp. The main 
results given in this paper do not depend on the occurence 
of the 'drag' effect, though components of the mass density 
matrix will appear in many of the equations. (The exception 
being that the 'drag' effect greatly increasing the scattering 
of electrons off the neutron vortices.) 



Throughout the rest of this paper all quantities will be 
understood to represent smooth averages over length scales 
large compared to the inter-vortex spacing but small com- 
pared to the length scales of physical interest (e.g., the stel- 
lar radius). In the smooth-averaged theory the effects of the 
underlying vortices manifest themselves as vortex forces act- 
ing on the smooth-averaged flow of the superfluids. Bekare- 
vich & Khalatnikov (1961) show a form for these forces can 
be found by allowing the energy density of the fluid to de- 
pend on the vortex densities. Mendell & Lindblom (1991) 
extended their method to mixtures of charged superfluids 
including the electromagnetic coupling. Following these ref- 
erences, first the vector fields A are defined so that the energy 
density due to the vortices, C/v, has a differential 



d(7v = A ■ d riv V- 

2m 



(13) 



It can then be shown that a vortex force per unit mass given 
by 



/v = -(V X A) X -— nvi^ 
p Am 



(14) 



acts on the smooth-averaged flow of the superfluid. 

It should be noted that the Bekarevich & Khalatnikov 
(1961) model ignores interactions between the vortices. Since 
the magnetic induction and circulation of a proton vortex 
falls off exponentially, ignoring the interactions between the 
proton vortices should be an acceptable approximation as 
long as the separation between these vortices is much larger 
than the London depth i.e., dp >> A. For simplicity inter- 
actions between neutron vortices will also be ignored here as 
including them would not significantly effect the results of 
this paper.* Finally, interactions between proton and neu- 
tron vortices, which can lead to pinning (or interpinning) of 
these vortices, must be considered. First note that Alpar, 
et al. (1984) have shown that the 'drag' effect causes the 
neutron vortices to carry a magnetic flux as well as the pro- 
ton vortices. The largest contribution to the pinning force is 
then due to the magnetic interaction of the vortices, giving 
a pinning energy per site of (e.g., see Mendell 1991a; Jones 
1991; Chau, Cheng & Ding 1992) 



1 

8^ 



ppp \ / 2$gcos6'„ 
ppp J \ 7rAsin6'„, 



10 ergs. 



(15) 



where 9np is the angle between the vortices. For d„ > dp the 
number of pinning sites per unit volume is roughly 



Afpin 



1 



d'^d 



dr, 



(16) 



The magnitude of the pinning force can be compared with 
the smaller of the vortex forces (that of the neutron vortices) 
by considering the ratio of the energy density associated 
with each force, i.e.. 



/pin ^ epinA''pin ^ „-4 j, 1/2 

— ~ ~ iU -O012 , 



(17) 



* Interactions between the neutron vortices can give rise to low- 
frequency oscillations of the vortex array known as Tkachcnko 
waves (Tkachenko 1966; Ruderman 1970; Fetter & Stauffcr 1970; 
Sonin 1987). 
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where e„ ~ 10^erg-cm~^ has been used (see MendeU 1991a). 
Thus pinning forces between the vortices will be ignored in 
this paper. 

Given the caveats of the last paragraph, expressions for 
An and \p are now found by constructing a model of the 
energy density. The magnetic energy density is comphcated 
by the fact that part of the magnetic induction is confined 
to the vortices, and also by the 'drag' effect which causes the 
neutron vortices to carry a magnetic flux, as noted above. 
Following Alpar et al. (1984) the magnetic induction due to 
the vortices is 



^oTi^pVp 



p„p 



'^on^nVn 



(18) 



Note that the first and second terms on the right side of this 
equation correspond to the magnetic induction due to the 
proton and neutron vortices, B^p and Bvn, respectively. 

In addition to the magnetic induction of the vortices a 
non-quantized magnetic field known as the London field, Bl, 
also exists whenever the vorticities of the fluids are nonzero 
(London 1960). Using equation (18), the London equations, 
(9) and (10), and B = Bvortcx + -Bl, gives 

B^ = -^I^u.p-^I^e^u:„. (19) 

I PV 

Ignoring interactions between the vortices, the vortex en- 
ergy density is just the vortex energy per unit length times 
the vortex density. The energy density of the vortices and 
electro-magnetic field then becomes 



u V ~\' Uem — SpTlvp -\- EnTivn 

1 2 

+ t; — ||Svortcx + Sl 
OTT 



Bt 



Bt 



D' 



(20) 



Note that the magnetic energy density of each vortex species 
is subtracted off since it has already been accounted for in 
the energy per unit length of the vortices [see eq. (4)] . It has 
also been assumed that the displacement field D equals the 
electric field E. Noting that electromagnetic energy density 
has the standard differential (e.g., see Jackson 1975) 



dt/E 



1 

47r' 



E dD 



4tv 



-H dB 



(21) 



and using equations (13), (18)-(21), and the London equa- 
tions, it can be shown that H = B, E = D, and 



An 



2m„ 



2m 






1 

An 



p u>p + apB 
h \uip + apB\ 



1 fmpC\'^ 
-[ — ) {u:p + apl 



(22) 



(23) 



Similar results for An and \p are given in MendeU (1991a). ^ 
The main purpose of this section is to present the 
equations for the neutron-star outer core given in MendeU 
(1991a) in the MHD limit, including the effects of the mag- 
netic field and the vortices (which are ignored in the final set 



' That paper includes the weak dependance of Sn on dn (which 
partly measures the interaction between these vortices) and sub- 
tracts i^vortex (which includes an interaction energy between the 
vortices) rather than just B^ + B^^ as in equation (20). Since 
interactions between vortices arc ignored in this paper, these dif- 
ferences are insignificant. 



of equations in that paper). Here the very small fraction of 
muons are ignored. Since the system is very degenerate the 
entropy density of the electrons and any excitations of the 
neutrons and protons are also ignored, while the largest of 
the dissipative terms found in MendeU (1991b) are included. 
The equations are for small perturbations away from a sta- 
tionary and axisymmetric equilibrium. (This allows weak 
interactions to be ignored.) The perturbed quantities are 
prefixed with a 5, and equilibrium quantities are given a 
subscript o. The equilibrium velocity Vo (vorticity u^o) is as- 
sumed to be that of uniform rotation with constant angular 
velocity while the equilibrium magnetic induction is taken 
to be an arbitrary stationary vector field Bo- 

The MHD equations then are as follows: the conserva- 
tion laws for the neutron and proton mass densities p„ and 
pp [note the appearance of the mass density matrix, given 
in eqs. (11) and (12)], given by 

dt5p„ + Vo ■ V5p„ + V ■ {p„„5v„ + p„p5vp) = 0, 



dtSpp + Vo ■ VSpp + V ■ {p„p5v„ + pppSvp) 



0, 



(24) 
(25) 



Maxwell's equations for the electric field E, magnetic in- 
duction B, charge density a, and current density J , given 
by 

epp epe 



0, 



Sa 



V ■ 5B = 0, 



V ■ (5B = 0, 



V X 5B = — dt5B, 

c 

Att 

V xSB = —5J, 

c 

o„ 5Ve ^ «o CD 

SE ^ X Bo X SB 

c c 

I — -5-Fn(mf) + <5-Fp{mf) ) , 

e \ Pp 



(26) 

(27) 
(28) 

(29) 
(30) 

(31) 



Newton's 2nd law for the neutron superfiuid velocity, v„, 
and the average velocity of the charged fluids, u, given by 

dtSvn = —Vo ■ V5v„ — 5v„ ■ VWo — V(5pn 



+ -^{Svp - 5v„) X iVo+S 
P„ 



-(V X An) X u:„ 



(32) 



dtSu — —Vo ■ V<5u — 3u ■ Vwo -'V3pp -'VSp^ (33a) 

Pp Pp 



+ — ( — X Bo ) + ^(5v„ - Svp) X ujo 

Pp \ C J Pp 



pp 



pp 



-(V X \p) X {tjp + apB) 



VS<p - P^5F^(^i) - —V • T, 

Pp Pp 



(33b) 
(33c) 
(33d) 



and Newton's law of gravitation for the gravitational poten- 
tial ip, given by 



WStp = AnGSp. 



(34) 



Preprint - Accepted for publication in MNRAS 



Mori. Not. R. Astron. See. 000, 000-000 (1997) Printed 1 February 2008 (MN plain T^ macros vl.6) 



Consider equations (26) and (27). These equations fol- 
low from the MHD assumption that charge neutrality is 
preserved, and equation (26) also relates p^ to Pp. Equa- 
tions (28) and (29) are just the usual Maxwell's equations, 
followed by equation (30) which is Ampere's law ignoring 
the displacement current term. This is justified for waves 
with frequency to and wave vector k when (cj/cfc)^ ^ 1, and 
this corresponds to the second part of the MHD limit. Next, 
equation (31) follows from considering the dynamics of the 
current density and is valid for frequencies small compared 
to the plasma frequency and in the limit of large conductiv- 
ity. Note that this equation implies that the Lorentz force 
on the electrons is vanishingly small in the MHD limit. Also 
in this equation, P„(mf) and Fpij^t) are the mutual friction 
forces due to electron scattering off the neutron and proton 
vortice respectively. These terms are small, though they are 
the largest of the dissipative corrections to this equation. 
Dissipative effects are discussed in Section 5. 

Now consider equations (32) and (33). In these equa- 
tions the /i's are the chemical potentials defined in equation 
(33) in Mendell (1991a). Equation (32) is just the usual Lan- 
dau equation for a superfiuid modified by the 'drag' effect 
and the vortex force of the neutron vortices [see eqs. (9) 
and (14)]. Equation (33) is for the average velocity of the 
charged ffuid defined by 



Pv , Pe 

Pp Pp 



(35) 



where pp = pp + p^, and it is the standard MHD equation 
for this quantity modified by the 'drag effect' and the vortex 
force of the proton vortices [see eqs. (10) and (14)]. Note 
that in equation (33) terms proportional to the equilibrium 
charge density, which is exceedingly small and in fact zero in 
non-rotating neutron stars, have been ignored [see eqs. (39) 
and (40) in Mendell 1991a]. The viscous stress tensor, r, and 
the mutual friction forces that appear in these equations are 
considered further in Section 5. 

One great simplification of the MHD limit is that the 
vector fields J, E, Vp, and Ve, are easily expressed in terms 
of B, u, and v„. Equations (30) and (31) already specify J 
and E and it can be shown that the velocities of the electrons 
and protons are given by 



SVe 



-Su - 



—Sv„ 



^SJ, 



epp 



epp 



(36) 



(37) 



pp pp 

ignoring factors of p^/pp in these equations only for the mo- 
ment. 

To finish this section it will be useful to review several of 
the characteristic frequencies associated with charged parti- 
cles and electromagnetic fields in terms of the angular fre- 
quency to, the wavevector k = (27r/A), and the wavelength 
A associated with plane wave solutions. 

Two frequencies well known in plasma physics are the 
plasma frequency, ojp, and the cyclotron frequency, ilc- For 
motion of the protons these frequencies are given by 



LUr = I 2 — ) ~ 3.2 X 10 ppi3 ' s 

fie = — ^ 9.6 X lO'^^Boia s-\ 



(38) 



(39) 



A characteristic frequency associated with charge currents 
in metals and superconductors is the helicon frequency, 
ojh, given by (see de Gennes & Matricon 1964; Abrikosov, 
Kemoklidze & Khalatnikov 1965) 



OJh 



fir 



C2fc2 



3.3 X 10 ppi3 As Boll s , 



(40) 



where Ae is A divided by 10*" cm. 

Of main interest in this paper are the Alfven waves 
found in the MHD limit of ordinary plasma physics. For an 
ordinary proton-electron plasma these waves travel along 
magnetic field lines with a frequency given by 



^Kp — V ficl^h ~ 



Bl/^-K 



Bo 



yiirp^ 



(41) 



These waves can be compared with those of a vibrating 
string with the tension of the field lines provided by mag- 
netic energy density (~ Bo/4jt) and the inertia provided by 
the mass that moves with the lines {pp). In the case of an 
ordinary-fluid neutron star the neutrons and protons would 
act as a single fluid and the density pp would be replaced 
with the total mass density p, giving 

Bo 



IjJa = 



=fc ^ 0.056pi,r^^^A6"^Boi2 s 



(42) 



where pis is the density divided by 10^^ g • cm^^. It will 
shown in the next section that Alfven waves do not occur in 
the exotic superfluid plasma considered here. 

There are also several characteristic frequencies associ- 
ated with oscillations of the underlying vortices. (See foot- 
note *; Krusius, et al. 1993; and Sonin 1987 for studies and 
reviews of vortex oscillations in superfluids; see Friedel, de 
Gennes & Matricon 1963; Tsui, et al. 1994; Kopnin, et al. 
1995; Blatter & Ivlev 1995; and Sonin 1996a,b for studies 
of vortex oscillations in superconductors.) The frequency of 
circularly polarized vortex oscillations which travel parallel 
to the vortices, known as Kelvin waves, are given by 



2mn Snk r,n„in-14 -1 \ -2 -1 

(^Vn ^ ^ ?5i 2.0 X 10 e„9Pnl5 K S , 



iOvp — 



h pp 



2.0 X 10 '^ep6Ppi3 'As 



(43) 



(44) 



for the neutron and proton vortices, respectively. In these 
equations p„i5 is p„ divided by 10^^ g • cm^'*, and £„9 and 
£p6 are £n and £p divided by 10^ erg-cm^^ and 10® erg-cm~^ 
respectively. These vortex frequencies can be shown to arise 
from the flrst terms in equations (22) and (23) for A„ and 



3 ALFVEN AND CYCLOTRON- VORTEX 
WAVES 

In this section is it shown that Alfven waves, which occur in 
ordinary plasmas, do not occur in the exotic plasma consid- 
ered in this paper. Instead a new frequency is shown to exist 
that depends on the cyclotron and vortex frequencies of the 
fluid. Thus it is referred to as the cyclotron- vortex frequency. 
However these waves are analogous to Alfven waves with the 
tension due to the magnetic energy density replaced by the 
vortex energy density. A straightforward derivation of this 
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frequency is given below. 

To begin with, consider uniform neutron-star matter. 
For simplicity assume that the matter is not rotating and 
that there is a uniform equilibrium magnetic induction 
which defines the z-axis, i.e., take = and B = BoZ. 
Next, look for plane wave solutions with the wave vector 
parallel to the z-axis, i.e., let k — kz. The space-time de- 
pendence of all the perturbed quantities can then be taken 
to be completely specified by exp(ifc2: — icut). Making these 
substitutions the equations in Section 2 become purely al- 
gebraic. The equations also decouple into components that 
are parallel and transverse to the z-axis. Ignoring dissipative 
effects, which will be considered in Section 5, the resulting 
equations for the transverse components are 



ik{z X 5E) = —SB, 
c 



An 
ik(z X SB) = —SJ, 
c 



SE = -(^ X Boz) 
c 



SVn = 0, 



(45) 



(46) 



(47) 



(48) 



ik 1 SJ 

—icoSu — — [{z X SXp) X apBoz] -\ ( — x Boz). (49) 

Pp Pp c 

A few steps will now show that Alfven waves do not 
exist in the solution to these equations. Using equation (46) 
it is easy to show that 



SB= ^(SJ xz) 
ikc 



(50) 



Then consider the vortex force acting on Su. Equation (23) 
implies that SXp is given by 



5Ai, 



2 2 

nipC 



h D,c 47re2 



[ik{z X Svp) + GpSB]. 



(51) 



Using equation (50) the vortex force in equation (49) can be 
simplified to become 



-^[(i; X 5\p) X apBoz] — -^{cuvp — u>h){Svp x z) 

Pp Pp 

I ujvp u)i ,sj 1 SJ 



(52) 



Note that the last term on the right side of equation (52), 
which arises from the interaction between the vortex and 
London fields, will cancel the last term in equation (49) 
[same as the first term in eq. (33b)]. Since this is the term 
that usually gives rise to Alfven waves these waves do not 
occur in the superconducting plasma considered here. It 
also should be clear that this result does not depend on 
the specifics of this plasma but should be true for any 
plasma with one component superconducting and in the vor- 
tex state. Furthermore this result applies not just to axial 
modes. In general the last term in equation (51) will give a 
force equal to (— l/47rp'p)(V x SB) x Bo which will give the 
same cancellation. 

From here on terms smaller by a factor p^/Pp will be 
dropped in all equations. Equation (36) already uses this ap- 



proximation, and substituting this into equation (47) gives 

(53) 



5E = -^ f ^ X BoZ 



pp \ c 

while equations (45) and (46) give 
4TTiuj 

From these one can show 

■^ iuo SJ 



!IhL(^^J- X Boz) 
epp c 



ppp{Su X z) 



U}„ 



c^k^ nl 



Bo + -p-{ — X Boz) 
c \lc c 



(54) 



(55) 



Substituting equations (37), (52), and (55) into equa- 
tion (49) and simplifying results in the following eigenvalue 
equation: 



z / \ rpp (-) 

^ — ^\\\yl^p — <^hj 5'c'^Vp 

Pp 



{SJ X Z) 



(56) 



-I- iuj{uJvp — 2u;h)SJ = 0. 
Solving equation (56), the eigenvalues are given by 

LU — LJhiuJvp — l^h) —0,cl^Vp = ±Uj{ljJvp — 2ct;h), (57) 

Pp 

and the components of the eigenvectors are related by 
SJx = ^^iSJy. The two sign choices in equation (57) give 
eigenfrequencies 



Ppp „ I^vp 

\lc(^vp r <^h 

Pp 2 



ppp o I "^VP 

ILcl^Vp H i^h 

Pp 2 



(58) 



where terms smaller than those kept by at least a factor of 
uJvp/Qc have been dropped. As far as this author knows the 
first term in equation (58) has not been derived previously 
from MHD theory. Since it depends on the cyclotron fre- 
quency fie and the Kelvin vortex oscillation frequency ci;vp 
it is natural to call it the cyclotron- vortex frequency, u>cv. 
In a neutron star it will have a magnitude given by 



(59) 




^{l + A)Boi2 



Note that ujcv is quite different from the Alfven frequency 
since it is nearly independent of the mass density (except 
through A) and depends on \/Bo- However, consider the def- 
inition of the critical field Hd = (<I>o/47rA'^) ln(A/^), above 
which the vortex state is favored in a type II superconduc- 
tor (and approximately the field within the London depth 
of one proton vortex). Noting that Sp — $oiifci/47r and 
nvp = Bo/^o then 



LJcv 



£pTl\rp 



k = 



HclBo/'ilT 



(60) 



It can seen then that ljcv is analogous to the Alfven fre- 
quency with the tension due to the magnetic energy density 
replaced by the vortex energy density SpUvp. This equals 
HciBo/4-K (roughly the magnetic energy density of the vor- 
tex array). Easson & Pethick (1977) found this same factor 



Preprint - Accepted for publication in MNRAS 



Mon. Not. R. Astron. Soc. 000, 000-000 (1997) Printed 1 February 2008 (MN plain T^ macros vl.6) 



in their study of the stress tensor for a type II supercon- 
ductor, though they did not consider its implication on the 
dynamics. 



4 INCLUDING ROTATION 

To see how the modes in the rotating case are affected by the 
cyclotron-vortex frequency, the solutions of the last section 
are extended here to the case of uniform rotation. Let the 
rotation axis define the z-axis of the system. To find plane 
wave solutions first transform to coordinates which corotate 
with the fiuid, given by 



X = X cos{Qt) + y sin{Qt) , 
y — ~xsin{Qt) + ycos{i^t), 



(61) 
(62) 



Substituting into equation (36) gives 



5J + ^(5j X i) = £^ P^Su + P^5v„ 
ujh trip y pp pp 

For uj ^ cjh the following approximations can be made: 

lOJ nip \ Pp Pp 

and from equations (37) and (50) 
Svp ^ 5u '-^Sv„, 



(73) 



(74) 



(75) 



a^SB^-(thh)(PlLSu+^5v„]. (76) 

\ U} J \Pp Pp J 

Substituting these into equations (68) - (71) and simplifying 



r =t. 



(63) 
(64) 



The spatial derivative operator V is covariant under this 
transformation, whereas the convective derivative of any 
fluid velocity V transforms as 

dtV + v„-W + V ■ Vvo = dt'V - 2n{V' x z). (65) 

If we also define new electric and magnetic fields 

(66) 



Vo 



E'^E+^xB 
c 



B' = B, 



(67) 



and ignore corrections of order {vo/c){uj/ck) <JC 1 then 
Maxwell's equations remain covariant under these transfor- 
mations. 

In accordance with the last section we consider axial 
modes with spacetime dependence &qy[ikz' ~ icut') and for 
simplicity assume the equilibrium magnetic field is aligned 
with the rotation axis. After making the above transforma- 
tions and dropping all primes one finds that Maxwell's equa- 
tions remain identical to those in the last section, and that 
the equations for the fluid velocities become 

^7^<^v„ - (Svn X z) ^ -^{Svp - (5n„) X z -\ (5A„, (68) 



—rrSu — {3u X z) = —^{5Vn — SVp) X Z 

212 pp 

+ 2npp ^^-+2?v:(t"^°")' 



where 



S\„ 



S\p = 



2Qh 



2 2 2 

nipC 
47re2 



(ikz X Sv„ 



(69) 



(70) 



2 2 

rUpC 



h{2Q + i}c) 47re2 



[ik{z X Svp) + UpSB]. (71) 



The goal will be to express everything in terms of 5v„ 
and Su. Using equations (31) and (54) one can show 



' epp J \LUi,J 



(72) 



icuSvn + 2Q^—^(5u X z) 

-h 20 ( h — :; 2 if (5t;„ X 2; = 0, 

2il CJp pnpi 



(77) 



20^— ^(5t;„ X z) + ^— ^ 1 Sv„ 

pp Pp \ LU LU 

+ 2n [ E^ + '^ - ^ - ^ ] (Su X z) 



+ iio 



Pp ' 2D. uj'^ 2D. 
iDcuJvp ppp 2iDcjh p 



(78) 



+ 

CO pp iO pp 



Su^O. 



To solve these equations one can invert equation (77) 
to solve for Su and then substitute into equation (78) to get 
a simple eigenvalue equation for Svn. For uj , D ^ uJh, <j-'vp, 
LUvn, and c^k^ ^ ojp the resulting dispersion relationship 
reduces down to 



uj-^±2D 



Pr. pp 

2 



Ppp \ 2 

LO 



+ ( 4n^-^ DcuJ^p— ) UJ T 2DDcuJ^p-^- = 0. 



Pr^Pp 



pp 



(79) 



p„pp 



where q^ = p„„Ppp — p'^p. The eigenvectors are related by 
Sux = liziSuy and 

Su=-^ (±^ ^) 5v„. (80) 

pnp y 2S2 pii / 

In the limit Bo = equation (79) can be solved exactly to 
give 



UJ = ± 



2D 



20- 



(81) 



Pr^pp 



where the upper branch is usually called the inertial mode 
and the lower branch is a new mode that appears due to 
the 'drag' effect. Keeping _Bo 7^ an approximate solution 
to equation (79) can be found by noting that pp <g; p„, 
p„„ « p„, and Q'^ « pnppp- The solution then is 

'20 



pp y pp V pp 



(82) 
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Note that the cyclotron-vortex frequency will dominate the 
effects of rotation only when 



B. > (pA (Hh^) (H.) 

\pp J \ e / \u;vpj 



2.4 X 10^ 
1 + A 



-ii2 '^6 VJ. 



(83) 



Thus this frequency will dominate the effects of rotation 
only in slowly rotating neutron stars with strong interior 
magnetic fields. In the more usual case we can expand the 
lower branch in equation (82) in powers of iUcv/'2^- Keeping 
only the lowest order terms, the result is 



UJ = ± 



'20 
20^ 

^CV Pp 

251 ppp 



(84) 



The top and middle branch are as in equation (81), while the 
bottom branch is a new low frequency mode that depends on 
iOcv Finally, recall that the frequencies in these equations 
are for the rotating frame; one must add =Ff2 to get the 
frequencies in the inertial frame. 



5 SIMPLIFIED EQUATIONS AND 
DISSIPATIVE EFFECTS 

It has been shown in the last two sections that of the vor- 
tex and magnetic forces the dominant term gives rise to 
cyclotron-vortex waves that replace Alfven waves. The next 
largest terms (e.g., consider ujvp/cocv) are roughly 15 or- 
ders of magnitude smaller for typical neutron star numbers. 
This suggests dropping the vortex force of the neutrons and 
approximating 5\p as 



SXr, 



2nip 



5B 



(85) 



It also will be convenient to assume that terms proportional 
to V(£p/|Bo|) are small and can be dropped. Terms smaller 
by a factor of p^j pp than those kept will also be ignored. 
Furthermore, equations (29) and (31) can be combined to 
eliminate &E from the equations, and equations (36) and 
(37) can be used to eliminate (5ue and &Vp. Finally, the small 
corrections to the velocity fields proportional to &J will also 
be dropped. With these approximations, equations (29) - 
(33) can be replaced by the following greatly simplified set 
of equations: 



dt5vn = — Wo ■ V5w„ — 5w„ • Vwo — V((5^„ + 5<^) 
-I- -^(l)u - Sv„) X uJo + 5F„(mf), 



(86) 



dtSB = V X 



(BiLSu + ^(5w„) xBo + VoXSB 

Pp Pp 



H — -l5f „(mf) + SFpl^f) 

ap \Pp 



(88) 



These equations are in a useful form for studing the effects 
of the interior magnetic field on the core spin-up during a 
pulsar glitch. More will be said about this in Section 6. 

The remainder of this section will discuss the dissipa- 
tive terms in these simplified equations. As discussed in 
Mendell (1991b) mutual friction due to electron-vortex scat- 
tering and the electron shear viscosity are the dominant dis- 
sipative effects (see also Cutler & Lindblom 1987; Cutler, 
Lindblom & Splinter 1990; Ipser & Lindblom 1991; Lind- 
blom & Mendell 1995). The viscous shear tensor is given 
by 5t — —rjeS®, where rje is the electron shear viscosity 
approximately given by 



V^ 



= 6.0 X lO^Vis^r.T^' g • cm" 



(89) 



where Ter is the temperature of the electrons divided by 
10^ K and S@ is the shear tensor, given by 



50% = i(V"5«e6 + VbSv: ^ iS^VJv^, 



(90) 



where a, b, c are space indices. The mutual friction forces are 
given in equations (14), (15), and (19) in MendeU (1991b). 
Keeping only the largest terms in those equations, and mak- 
ing the same approximations as discussed at the beginning 
of this section, one can show that 



5-F'„(mf) — — B,^ 



5f p(mf) — ~Bp 



P^Pp 



kol <^ 5V 



CJo(i^o ■ SV) 1 

J 



h^o 



pp<i>o 



(V X SB) 



(91) 



(92) 



where SV = 5v„ — Su. A term proportional to Bo ■ (V x 
SB) has been dropped in equation (92) in keeping with the 
approximate form for SXp given in equation (85). Simple 
estimates for the mutual friction coefficients based on the 
formulas given in Mendell (1991b) are 

B„ ^ ■Ol{p,^p/ ppp)^ {ppp/ Pp)'^^^ {pp/ Pp.) 



10' 



B„ 



M{pjp,) ^ 10- 



(93) 
(94) 



It will now be possible to find an energy functional from 
which the damping time r of modes with time dependence 
given by exp{iujt — t/r) can more easily be found. General- 
izing equation (32) given in Mendell (1991b) an appropriate 
real energy functional is 



B = - / < p„„5u„ • SVn + PppSVn ■ Su* 



dfSu — — Vo ■ 'VSu — Su ■ Vuo — ^{Sfip H -Sfj,e 

nip 

— VS(p + ^—^{SVn — Su) X LJo 
pp 

+ /^„ , (V xSB) X Bo 



Pp^o\Bo\ 



Pr. 



■SF 



1 



nfmf) ^ ■ T, 

pp 



(87) 



-I- pnpSu ■ Sv^ + pppSu ■ Su* + y 



_dPUi_ 

dpAdpB 



SpASp*B (95) 



VStp ■ VStp* EpSB ■ SB* 



47rG 



^o\Bo\ 



xTx. 



In this equation A and B are indices which take on the 
values (n,p, e) and Ui is the internal energy of the fluid 
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{dUi/dpA = ma)- The damping time r of a mode can then 
be calculated using 



(96) 



1 ]_(1E 

7 ^ ~2E~dt' 

Using the equations in Section 2 and the approximations of 
this section to compute dE/dt, and writing the result in the 
form 1/r — 1/r^ + 1/rmf, the shear damping time and the 
mutual friction damping time can be shown to be 



1 



1 
2E 



2r]e5Q°' SQaijd' X, 




1 I 4 



\sv[ 



\u^o-SV\- 



(97) 



(98) 



B 



pp$Sfflp|-Bo| 



The damping times of the modes found in Section 3 and 
Section 4 can now be estimated. First, equations (80), (36), 
and (76) are used to relate 5vn, Sve, and SB to 5u, and the 
results then are substituted into equations (97) and (98) . It 
will be convenient to express the damping times in terms of 
the quality of the mode, Q — u)t. For the cyclotron-vortex 
frequency, given in equation (59), the viscous damping time 
dominates, and Qcv = 0Jcv{pp/ ppp)[pp/{j)ek'^)], or 



Ocv ~ 2.7 X lO^AeT.r^is"' f ^ j ^(1 + A)Boi2 . (99) 

For the rotating case, the frequency splits into the three 
branches given in equation (84), and viscosity dominates the 
top branch, while mutual friction due to electron scattering 
off the neutron vortices dominates the damping of the mid- 
dle and bottom branches. The resulting Q for each of these 
modes is 

' ^P ^ in4o -1\ 2Tn 2 
-2 « 10 ihpis Ae le7 

^ = \ [B„(p„/p,)]-^ ^ WO{p„/p^pf{pJp„)'/^ (1°°) 

.[B„(p„/p,)]-i ^ 100(Pp,/p„,)2(Pp/Pp,)1/2 



6 MAGNETIC SPIN-UP OF THE CORE 

It has been suggested that magnetic fields are responsible for 
the spin-up of the neutron star core during a pulsar glitch. 
Detailed studies of the spin-up process are given in Easson 
(1979), Chao, Cheng & Ding (1992), Sedrakian, et al. (1995) 
and Abney & Epstein (1996). It also has been suggested that 
Alfven waves might be excited during the glitch (cf., Baym 
et al. 1969; Van Horn 1980 and references therein.) The goal 
of this section is to study the influence of cyclotron-vortex 
waves on the crust-core coupling time, and to place limits 
on the interior magnetic field. First the simple method of 
Abney et al. (1996) is used, followed by estimates based on 
the method of Easson (1979). The results are compared with 
the other studies mentioned above. 

Recently Abney et al. (1996) have shown that the rapid 
spin-up time for the 1988 'Christmas' glitch in the Vela pul- 
sar implies a crust-core coupling time of less than 10 sec- 
onds. They also argue, following Abney & Epstein (1996), 
that Ekman pumping and viscosity cannot be responsible for 



this short timescale. Following their 'back of the envelope' 
derivation, any wave with phase velocity v — uj/k associ- 
ated with the forces that spin up the core must obey the 
approximate constraint 

R Rk 



<tci, 



(101) 



where fci is the crust-interior coupling time and R the radius 
of the core. This equation can be used to derive a lower limit 
for the interior magnetic field. 

For an ordinary-fluid neutron star, substituting the 
Alfven wave frequency given in equation (42) into equation 
(101) and solving for Bo gives a lower bound on the interior 
magnetic field of 



Bo > lO'Vis'^'-Reiciio"' G, 



(102) 



where Re is R divided by 10^ cm and tcno is tci divided by 
10 s. This is exactly the same lower bound that Abney et 
al. (1996) give for an ordinary-fiuid neutron star. 

However it has been shown here that Alfven waves do 
not exist in superconducting-superfiuid neutron stars, but 
are instead replaced by cyclotron- vortex waves. Substituting 
equation (59) into equation (101) and solving for Bo gives 



Bo > 



W 



1 + A 



Pp \ n 2, -2 p 



(103) 



To compare with Abney et al. (1996), in this case they 
assume that confining the flux to the proton vor tices in- 
creases the Alfven phase velocity by a factor of y^Hci/Bo 
[equals the square root of the vortex energy density over 
the magnetic energy density, see eq. (60) and Easson & 
Pethick 1977] and use Hoi ~ 10^^ G. This gives roughly the 
cyclotron-vortex wave velocity. Thus their result is funda- 
mentally the same as that given here, though they have ig- 
nored the mass density dependance of Hci and do not flnd 
the factor (pp/ppp)'^- 

To improve on the simple method of Abney et al. (1996) 
one should account for the effects of rotation, compress- 
ibility, stratification, dissipation, gravity, and geometry on 
the spin-up. Previous studies have shown that the spin-up 
of a rotating fluid proceeds by boundary-layer formation, 
followed by suction into the boundary layer (e.g., Ekman 
pumping, as mentioned above), and finally damping of any 
residual oscillations (Easson 1979; Abney & Epstein 1996; 
and references therein). Easson (1979) considered the spin- 
up of a constant-density rotating-slab model using MHD 
equations for an ordinary p-e plasma, and concluded that 
suction into a boundary layer was the most important pro- 
cess, even in the case of magnetic spin-up. The result given 
in that paper for the magnetic spin-up time of the core de- 
pends on both the angular velocity and the Alfven frequency, 
and is given by {QkR/ujAf^'H'^/^)- However, Abney & Ep- 
stein (1996) have shown that stratification (and compress- 
ibility) in a neutron star can limit the spin-up during the 
suction phase to only a very small region near the crust- 
core boundary, greatly increasing the spin-up time of the 
rest of the interior fluid. That paper only considered the ef- 
fects of viscosity on the spin-up, but its general conclusions 
should apply in the magnetic case as well. Easson's mag- 
netic spin-up time then would only apply to the crust-core 
boundary-layer, which he found has an approximate thick- 
ness oiuiA/i^k) (the distance an Alfven wave penetrates the 
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core during one rotation). For the entire core to spin- up, one 
can then imagine a succession of boundary layers forming 
from the crust down to the centre of the star. Each layer 
would spin up on a time scale approximately equal to Eas- 
son's magnetic spin-up time and the number of such layers 
would be QkR/uA- This gives a total magnetic spin- up time 
of (QkR/cdA)^ (2/^). Obviously this represents an upper 
bound on the spin-up time, and Easson's result represents a 
lower bound. To summarize, if magnetic fields are responsi- 
ble for the core spin-up, one can tentatively predict that for 
an ordinary-fluid neutron star the crust-core coupling time 
will be constrained by. 



V Wa / fl ~ "^^ ~ \ uJa J ^ 



(104) 



Ignoring dissipative effects and the dynamics of the neu- 
trons, it can be shown that the simplified MHD equations 
given at the beginning of Section 5 are identical to those 
in Easson (1979) but with the magnetic force in that paper 
increased by (wcv/wa)^ and the velocity field in that pa- 
per replaced with {ppp/pp)u. Thus the corresponding limits 
on tci in the superconducting-superfluid case are identical 
to those in equation (104), but with the Alfven frequency 
replaced with the cyclotron- vortex frequency, i.e., 

-flkRy'^'^ 2 



^"^^Y'^^<.c:<(^)"''l. (105) 



V Wcv 

Easson correctly predicts the same lower limit given here (up 
to factors of Ppp/pp) by assuming that the Alfven phase ve- 
locity is increased by a factor of y^Hci/Bo (the same factor 
used by Abney et al. 1996 as described above). 

Using equation (42) for the Alfven frequency and equa- 
tion (59) for the cyclotron-vortex frequency in equations 
(104) and (105) gives the following limits on the interior 
magnetic field: 



10 



12 Pv. 



1/2 



Re 



O 1/2+ 3/2 

in the ordinary fluid case, and 



^ 1/2 p o 2/5 



.3/5 



(106) 



10** pI 



Re. 



1 + A p2^ ^2*0110 



G<B<, < 



10^ 



^2 p2(-> 4/5 

Pp Re^2 (^Q7) 



l + Aplp tcno6/5 



in the superconducting-superfluid case. Note that the lower 
bounds in these equations are much less stringent that those 
given by Abney et al. (1996), and also in equations (102) and 
(103) (e.g., consider the Vela pulsar with Q, — 70.6 s~^ and 
tci close to 10 s). But if stratiflcation is as important in the 
case of magnetic spin-up as the study of Abney & Epstein 
(1996) implies, then the lower bounds of Abney et al. are 
likely to be correct. A more rigorous study of the spin-up 
process is needed to confirm the upper bounds given above, 
and may in fact show that the lower bounds should actually 
be very close to these upper bounds. 

Without direct observation it is hard to say much with 
certainty about the interior magnetic field of a neutron star. 
An early estimate predicted it would persist for the life- 
time of the universe (Baym, Pethick & Pines 1969), but 
other studies suggest that it might decay in less than lO^yrs 
(Jones 1987; Haensel, Urpin & Yakovlev 1990) or that pin- 
ning between neutron and proton vortices causes a portion 
of the interior flux to be expelled during spin down (Srini- 
vasan et al. 1990; Jones 1991; Ding, Cheng & Chau 1993). 



These latter models of flux expulsion give results that are 
consistent with the lower bound given in equation (107), but 
arc harder to reconcile with the lower bounds given in the 
other equations. 

In contrast to the lower bounds given here, Chao, Cheng 
& Ding (1992) actually find an upper bound on the interior 
magnetic field of approximately Bo < 10*~^G using vortex 
creep theory to study the core spin-up. They assume that 
the neutron vortices are pinned by their magnetic interac- 
tion with the proton vortices, but as they point out [and see 
eq. (17)] other forces dominate the pinning forces during the 
initial spin-up. It is possible that pinning and vortex creep 
could be important during the transition from spin-up to 
post-glitch relaxation back to equilibrium. In even further 
contrast to the results given here, Sedrakian et al.(1995) use 
a vortex cluster model which predicts that proton vortices 
cluster around the magnetized neutron vortices. This pro- 
duces an interior magnetic field that is tied to the star's 
rotation rate and in their model the spin-up is independent 
of any relic field present at the star's birth. 

It is hoped future observations of pulsar glitches that 
time-resolve the spin-up will help distinguish between the 
various models mention here. 



7 DISCUSSION 

The MHD equations for the outer-core of a superconducting- 
superfluid neutron star have been presented in equations 
(24)-(34), while a greatly simplified form of the equations 
for the vector fields is given in equations (86)-(88). It has 
been shown that Alfven waves are replaced by cyclotron- 
vortex waves in this region, and it has been shown how these 
waves and the 'drag' effect modify the inertial mode when 
rotation is included [see eqs. (59) and (84)]. The dissipative 
time scales due to electron viscosity and mutual friction can 
be computed using equations (95)-(98), and the damping 
times of the modes discussed in this paper can be found 
in equations (99)-(100). The magnetic spin-up of the core 
during a pulsar glitch is discussed in Section 6, and various 
limits on the interior magnetic field are found and compared 
with other results. 

It is intriging to think that glitch observations can con- 
strain the physics of the neutron-star core. The limits given 
in Section 6 should be improved upon by solving the MHD 
equations given in this paper for the core spin-up, either by 
using the analytic methods of Easson (1979) or Abney & Ep- 
stein (1996), or numerically. The effects of the magnetic field, 
rotation, compressibility, stratification, dissipation, gravity, 
and geometry on the spin-up could then be studied. Other 
areas that require further investigation are the crust-core 
boundary conditions (including surface pinning effects) , the 
effects of having an inner core, and the relativistic correc- 
tions to the equations. The results of future studies could 
then be compared with the competing models mentioned 
in Section 6, and ideally it will also be possible to make a 
comparison with future observations. 
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